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Abstract 

The short-time asymptotic behavior of option prices for a variety of models with jumps has received much 
attention in recent years. In the present work, a novel second-order approximation for ATM option prices under an 
exponential tempered stable model, a rich class of Levy processes with desirable features for financial modeling, is 
derived and, then, extended to a model with an additional independent Brownian component. Our method of proof 
is based on an integral representation of the option price involving the tail probability of the log-return process 
under the share measure and a suitable change of probability measure under which the process becomes stable. 
Our approach is sufficiently general to cover a wide class of Levy processes which satisfy the latter property and 
whose Levy densities can be "closely" approximated by a stable density near the origin. The results hereafter shed 
new light on the connection between both the volatility of the continuous component and the jump parameters 
and the behavior of ATM option prices near expiration. In the presence of an additional Brownian component, the 
second-order term, in time-t, is of the form di i' 3 ~ y ^ 2 , with the coefficient <fe depending only on the overall jump 
intensity of the process and the tail- heaviness parameter Y. This extends the known result that the leading term is 
at 1 / 2 /y/2n, where a is the volatility of the continuous component. In contrast, under a pure-jump tempered stable 
model, the dependence on the overall jump intensity and Y is already reflected in the leading term, which is of the 
form dit 1 ' Y . The information on the relative frequency of negative and positive jumps appears only in the second- 
order term, which is shown to be of the form dit and whose order of decay turns out to be independent of Y. The 
asymptotic behavior of the corresponding Black-Scholes implied volatilities is also addressed. Our numerical results 
show that first-order term typically exhibits rather poor performance and that the second-order term significantly 
improves the approximation's accuracy. 
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1 Introduction 

It is generally recognized that the standard Black-Scholes option pricing model is inconsistent with options data, while 
remaining a widely used model in practice because of its simplicity. Exponential Levy models generalize the classical 
Black-Scholes setup by allowing jumps in stock prices while preserving the independence and stationarity of returns. 
There are several reasons for introducing jumps in financial modeling. First of all, sudden sharp shifts in the price level 
of financial assets often occur in practice, and this "jump" risk is quite hard to handled within continuous-paths models. 
Second, historical asset prices exhibit distributions with so-called stylized features, such as heavy tails, high kurtosis, 
volatility clustering and leverage effects, which are hard to replicate within purely-continuous frameworks. Finally, 
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market prices of vanilla options exhibit skewed implied volatilities (relative to changes in the strikes), in contrast to 
the classical B lack-Scholes model which predicts a flat implied volatility smile. Moreover, the fact that the implied 
volatiliy smile and skewness phenomenon becomes much more pronounced for short maturities is a clear indication of 
the presence of jumps. We refer the reader to Cont and Tankov [5] for further motivations on the use of jump processes 
in financial modeling. 

One of the first applications of jump processes in financial modeling is due to Mandelbrot (25], who suggested a pure- 
jump stable Levy process Z to model power- like tails and self-similar behavior in cotton price returns. Merton [31j and 
Press 34 subsequently considered option pricing and hedging problems under an exponential compound Poisson process 
with Gaussian jumps and an additive independent non-zero Brownian component. A similar exponential compound 
Poisson jump-diffusion model was more recently studied in Kou [55], where the jump sizes are distributed according to 
an asymmetric Laplace law. For infinite activity exponential Levy models, Barndorff-Niclsen [T] introduced the normal 
inverse Gaussian (NIG) model, while the extension to the generalized hyperbolic class was studied by Eberlein, Keller 
and Prause [10]. Madan and Seneta [28] introduced the symmetric variance gamma (VG) model while its asymmetric 
extension was later analyzed in Madan and Milne [27J and Madan, Carr and Chang [23]. Both models are built on 
Brownian subordination; the main difference being that the log-return process in the NIG model is an infinite variation 
process with stable like (a — 1) behavior of small jumps, while in the VG model the log-price is of finite variation with 
infinite but relatively low activity of small jumps. The class of "truncated stable" processes was first introduced by 
Koponen [24]. It was advocated in financial modeling by Cont, Bouchaud and Potters [7J and Matacz [30], and further 
developed in the context of asset price modeling by Carr, Geman, Madan and Yor [5j, who introduced the terminology 
CGMY. Nowadays, the CGMY model is considered to be a prototype of the general class of models with jumps and 
enjoys widespread applicability. 

The CGMY model is a particular case of the more general KoBoL class of Boyarchenko and Levendorksii [4], which 
in turn is a subclass of the semi parametric class of tempered stable processes (TSP) introduced by Rosiiiski [39]. TSP 
form a rich class with appealing features for financial modeling. One of their most interesting features lies in their 
short-time and long-time behaviors. Denoting by Y G (0, 2) the index of a TSP, L :— (L t )t>o, it follows that 

(^- 1 L, lt ) t>Q A (Z f ) t > , (h^O), (1.1) 



where Z := (Z t ) t>0 is a strictly a-stable Levy process with parameter a = Y. In fact, (1.1 1 holds when Y G (0, 1) 
provided L is driftless, while it holds when Y G (1, 2) regardless of the "center" or mean of L. Roughly, ( 1.1 1 suggests 
that the short-time behavior of TSP is akin to that of a stable process, with its heavy-tailed distribution and self- 
similarity properties which are desirable for financial modeling as emphasized, for example, by Mandelbrot |29j . In 
contrast, 

{h- x L ht ) t ^ A (W t ) t > , (fc->oo), 

where (Wt) t>0 is a Brownian Motion, which suggests that in the long-horizon, the process is Brownian-like. In terms of 
increments, if we were to consider the consecutive "high-frequency" increments of the process L, say \Lh% ~ ^h{i~\)\%>\ 
with h <C 1, these will exhibit statistical features consistent with those of a selfsimilar stable time series. But, low- 
frequency increments, say {L^ — Lh(i-i)}i>l with ft 3> h, will have Gaussian like distributions. 

Stemming in part from its importance for model calibration and testing, small-time asymptotics of option prices 
have received a lot of attention in recent years (see, e.g., [2]. [3]. [TT]. [12]. [15] . [T7]. [15] . [T9]. [20]. [2"T]. [33]. [35]. [43]). 
We shall review here only the studies most closely related to ours, focusing in particular on the at-the-money (ATM) 
case. Carr and Wu |jJj first analyzed, partially via heuristic arguments, the first order asymptotic behavior of an Ito 
semimartingale with jumps. Concretely, ATM option prices of pure-jump models of bounded variation decrease at the 
rate 0(t), while they are just Oit 1 ^ 2 ) under the presence of a Brownian component. By considering a stable pure-jump 
component, [9j also showed that, in general, the rate could be 0(t beta ), for some /3 G (0, 1). Muhle-Karbe and Nutz [52] 
formally showed that, under the presence of a continuous-time component, the leading term of ATM option prices is of 
order \/t, for a relatively general class of Ito models, while for a more general type of Ito processes with a-stable-likc 
small jumps, the leading term is Oil; 1 / 01 ) (see also [T5J Proposition 4.2], [TSJ Theorem 3.7], and [15] Proposition 5] for 
related results in exponential Levy models). However, none of the these papers have obtained second or higher order 
asymptotics for the ATM option prices, which are arguably more relevant for calibration purposes, given that the most 
liquid options are of this type. 
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In the present paper, we study the small-time behavior for at-the-money call (or equivalently, put) option prices 



E (S t - S Q y = S Q E [e Xt - 1 



(1.2) 



under the exponential Levy model 



S t ■= S e Xt , 



(1.3) 

where X := (X t )t>o is the superposition of a proper tempered stable Levy process L := (L t )t>o in the sense of |39j (see 
Section [2] below for a precise definition) and of an independent Brownian motion (o~W t )t>o] i- e -j 



X t :=(rW t +Lt, 



(1.4) 



where W :— (Wt)t>o is a standard Brownian motion independent of L. Throughout, x + is the positive part of x. As 
it will be shown below, the first order asymptotic behavior of (1.2) in short-time under (1.4) takes the form: 



limt-^ Y E(S t - S Q )+ = S Q E{Z+), 



i-+0 



(1.5) 



where Z is a centered a-stable random variable with a = Y under P. When a 7^ 0, Z ~ JV(Q,a 2 ) (i.e., a = 2) and, 
thus, E(Z + ) = a/^/2n, as already shown in [43] and [35]. When cr = and L is symmetric, the characteristic function 
of Z is explicitly given by 

Ee iuZ = e -2cr(-Y)\cos(iY7r)\\u\ Y ^ 



where C is a measure of the overall jump-intensity of the process. In that case, (see (25.6) in [42]), 



di := E(Z + ) 



2CT(-Y) 



ttY 



l/Y 



(1.6) 



We refer the reader to Remark |3 . 5 1 below for the explicit formula for the leading term in the general case. Interestingly 
enough, under the presence of a continuous component (a =/= 0), the first-order asymptotic term only reflects information 
on the continuous-time volatility, in sharp contrast with the pure-jump case where the leading term depends on the 
overall jump-intensity parameter C and the index Y , which in turn controls the tail-heaviness of the distributions. 



The asymptotic formula (1.5) is in agreement with the result of Tankov [43] (see Theorem 5 therein), who showed 



that for a pure-jump Levy process, whose characteristic function is similar to that of a stable process, the rate of 
convergence of ATM option prices is of order t x l Y . Concretely, it is shown in [?3] that if the characteristic function of 
the Levy process (X t )t>o is such that 



E e 



iuX± 



exp (iuj - \u\ Y f(u)) 



(1.7) 



for 1 < Y < 2 and a function / satisfying lim tt /'oo f( u ) = c+ and lim^v^-oo f{u) = c_ for some constants < c+, c_ < 
00, then 



hm<- 1 / Y E(^-S* ) + -S'o^ 



c_ 



(1.8) 



As it turns out, the characteristic function of a tempered stable process satisfies the specification (1.7) (see Theorem 
2.9 in [39 ) and, thus, (1.8) can be applied to explicitly write the leading order term in (1.5). Alternatively, we can also 
use a result of |32j (see Theorem 4.4 therein), where a formula similar to ( 1.5 ) is shown for a more general class of pure- 



jump martingales. Nevertheless, as previously mentioned, none of the these papers obtained second or higher order 
terms for the ATM option prices and the approaches therein do not seem amenable to study higher-order asymptotics. 



The main result of the present paper establishes a second order correction term for the approximation (1.5) in a 



general tempered stable Levy model with a possible nonzero independent Brownian component. To the best of our 
knowledge, these are the first results of this kind in the literature of exponential Levy models. We show that the 



second-order asymptotic behavior of the ATM call option price (1.2) in short-time is of the form 



1 

5o 



E(S t -So) + = dif +d 2 t + o{t), (t->0), 
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in the pure-jump tempered stable Levy case (i.e., a = 0), while in the presence of a non-zero independent Brownian 
component (i.e., a ^ 0), 

-^E(S t - S ) + = dit* + d 2 t^ + o (t 3 -^) , (t -> 0), 
bo v J 

for different constants d\ and that we will determine explicitly. To wit, we found that, under the presence of a 
nonzero Gaussian component, the second-order term depends only on the overall jump intensity parameter C and the 
tail-heaviness parameter Y . However, for a pure-jump case, the tail features of the Levy measure are already present 
in the second-order term. The above asymptotic behaviors should also be compared to the corresponding behavior 
under the standard Black-Scholes model, where it is known that (see, e.g., O Corollary 3.4]) 

E ( e ^-4* _ i)+ = * - -^=ti + o(th. 

V2^ 24 v / 2^ 

As a byproduct of our asymptotic results for option prices, we also give the asymptotic behavior of the corresponding 
Black-Scholes implied volatilities. 

The derivation of our second-order approximation results builds on two facts. First, as in 13], we make use of the 
following representation of Carr and Madan [B] : 

1 f°° 

—E(S t ~So) + =F*(X t >E)= e- x F*(X t >x)dx, (1.9) 

where P* is the martingale probability measure obtained when one takes the stock as the numeraire (i.e., ¥*(A) :— 
£(5^1^)) and E is an independent mean-one exponential random variable under P*. The measure P* is also called 



the share measure (see [B]). As it turns out, under P*, (Y t ) t > admits a decomposition similar to (1.4 1 



X t = crW t * +L* t , t>0, (1.10) 

where W* := (W*)t>o is a Wiener process while L* := (L*) t >Q is a proper tempered stable process, independent of 
W* . Second, we change probability measures from P* to another probability measure, say P, under which (££) t >o is 
a stable Levy process and (W t *)t>o is still a standard Brownian motion independent of L* . This change will in turn 
enable us to exploit some key features of stable processes such as self-similiarity and the tail behavior of their marginal 
densities. 

As explained above, the class of tempered stable Levy processes is quite rich and enjoys several desirable properties 
for financial modeling. A natural question is whether our approach is general enough to handle any Levy process, which 
satisfies the fundamental property of being stable under a suitable change of probability measure. As a byproduct of 



our analysis (see Remark 4.4 for precise arguments), we found this to be the case under a key additional condition 
implying that the Levy density of the process can be "closely" approximated by a stable density near the origin. 
Concretely, denoting the Levy density of the process under the risk- neutral measure by s, our results will hold true 
(modulus some additional mild conditions) if, for some Y £ (1,2), the function q(x) := s(x)/\x\ Y+1 is such that 

lim - {C+ - q(x)) = M, lim - (C_ - q(x)) = —G, (1.11) 

x\0 X i/O X 

for some constants C+, C_, G > and M > 1. 

The present article is organized as follows. Section [2] contains preliminary definitions and properties related to 
tempered stable processes and probability measure transformations, which will be needed throughout the paper. Section 
3 establishes the second-order asymptotics of the call option price under the pure-jump tempered stable model (a = 0). 
Section 4 establishes the second-order asymptotics of the call option price under the tempered stable model with an 
additional independent non-zero Brownian component (cr 7^ 0). Section [5] illustrates our second-order asymptotics for 
the important particular class of CGMY models, recovering our preliminary results first presented in [16|. Section [6] 
assesses the performance of the asymptotic expansions through a detailed numerical analysis for the CGMY model. 
The proofs of our main results are deferred to the Appendices. 
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2 Tempered stable processes 



Let L := (L t )t>o be a pure-jump Levy process with triplet (0, b, v) and let W :— (Wt)t>o be a Wiener process, 
independent of L, defined on a complete filtered probability space (f2, J 7 , (J 7 t)t>o> P)- We assume zero interest rate and 
that P is a martingale measure for the exponential Levy model 



S t = S e 



with X t := aW t + L t 



where (St)t>o represents the price process of a non-dividend paying risky asset. Equivalcntly, the Levy triplet (tr 2 , b, v) 
of X := (X t )t>o is such that 



1 =E e 



exp t 



bu+Y+l { eX - 1 - xl {\x\<i}) v{dx) 



(*>0), 



(2-1) 



where hereafter H,q := R\{0} and the Levy triplets are given relative to the truncation function l{|a;|<i} (see Section 
8 in [H]). Without loss of generality, we also assume that (X t )t>o is the canonical process X t (iS) := uj(t) defined 
on the canonical space fl = O([0, oo),R) (the space of cadlag functions ui : [0, oo) — > H.) equipped with the cr-field 
J- := <t(X s : s > 0) and the right-continuous filtration J- t := C\ s> ta{X u : u < s)). 

In this paper, we consider an exponential tempered stable Levy model, i.e., throughout we assume that, under P, 
L := (L t )t>o is a proper tempered stable Levy process in the sense of |39j with index Y G (1,2). That is, the Levy 
measure v of L under P admits a density s : Rq — ► [0, oo) of the form (see (2.7) in [39]): 



s(x) = \x\ Y 1 q{x), 



with a function q : R\{0} — > [0, oo) given by 

q(x) := 

for a nonzero measure p on Rq such that 



q + (x) x > 0, 
q~{— x) x < 0, 



/ °° e- x ' s s Y p(ds) x>0, 
f°e-*/°\s\ Y p(ds) x<0, 



(2.2) 



\s\ Y p(ds) < oo. 



By Bernstein's theorem, the representation ( 2.2 ) is equivalent to asking that q + : R + — > [0, oo) and q 



.+ — > [0, oo) are 

completely monotone functions, meaning that (— l) n j-^q ± (x) > 0, for all x > 0. In particular, q ± : R + — > [0, oo) arc 
strictly decreasing and convex. A prototypical tempered stable process is the CGMY process (also called the KoBol 
model) of widespread applicability in mathematical finance. For the CGMY process, p(ds) := CM Y 8{M- 1 }{ds) + 
CG Y 8s_q-i\ (ds) so that the Levy measure of L is given by 



v{dx) 



Ce~ Mx Ce Gx , , 

X l+Y 1 U>o} + yJT+Y X {^<o} I dx - 



(2.3) 



For P to be a martingale measure, we need that Ee^* < oo or, equivalently, that 



OO POO 

e x v(dx) = / e x q(x)x~ Y ~ 1 dx < oo. 
i Ji 



It is not hard to see that a necessary and sufficient condition for (2.4 1 to hold is that 

p((l,oo)) = 0. 



Indeed, suppose that (2.4 1 holds. Then, for any sq > 1, we will have that 



I e x q(x)x Y 1 dx > J 



-Y-l 



- x d-Vs Y p(ds)dx > J X 



OO p OO 



(2.4) 



(2.5) 



poo 

dx / s Y p(ds) = oo, 

J Sn 



■5 



unless p([sq, oo)) = 0, for any s > 1. On the other hand, if (2.5) is satisfied, q(x) = f e x / s s Y p(ds) < e x J Q s Y p(ds) 
and, thus, 

1 />oo />1 />00 

e x u(dx) = / e x q(x)x~ Y ~ 1 dx < / s Y p(ds) / x" Y_1 dx < oo. 



Therefore, from now on, we shall enforce the condition (2.5). We also make the following standing assumptions: 



Assumption 2.1. The function q : M. — >• [0, oo) is such that 

In q \Xj\ 

(i) limsup — q'(x) < oo, (ii) limsupg (x) < oo, (iii) limsup — j— ; — < oo. 

x\0 i/0 |a|->oo \ x \ 



(2.6) 



Since 



-q'{x) = 



,-x/sY-l 



s y - 1 p{ds) > 0, for x > 0, g'(x) = 



-x/a\ a \Y-l 



s\ Y - L p(ds)>0, {otx<0, 



conditions (2.6 i) and (2.6 -ii) are equivalent to asking that the measure p satisfies 

f \s\ Y - 1 p(ds) <oo. 



Hereafter, we set 



C + :=q(0 + ) = / s Y p(ds), C- :=g(Q-) 
Jo 



\s\ Y p(ds) 



It is also convenient to introduced the following standardized versions of q and p: 
p(dx) :■ 



1 1 

^l{x>0} + ^-l{a;<0} ) P(dx), 



-i \ l( x ) -, l( x ) i 

Q\ x ) l{x>o} + -g— 1 { K <o} 



J °° e- x ' s s Y p{ds) x > 0, 
/° e-*/»| a |^(ds) z<0. 



Note that these standardizations ensure that 



q(0+) = / S r p(ds) = 1, g(0") = / \s\ Y p(ds) = 1. 

JO J-oo 



(2.7) 



(2.8) 



(2.9) 



Under (2.9), hereafter q(0) := 1. The following properties of q are crucial to obtain the second-order asymptotics: 
Lemma 2.2. Let 



M := -q'(0 + ) = / s Y - 1 p(ds), and G := g'((T) 
Jo 



kl y - x p(d S ). 



Then, under the conditions \2.6\i) and (2.6 -ii) 



1 



Ivo 5 ^ ~ = yM1 {v>°} + \y\ G1 iv<o}- 



(2.10) 



In particular, 



(i) Mm ?(fy) 1/4 = e-v M l {v>0} + e-^ G l {v<0} , (ii) lim ] lnq(5y) = -yMl {y>0} - \y\Gl {y<0} . (2.11) 
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Proof. Let y > 0. Then, for 5 > 0, 

1 



(i-g(fy)) 



where the limit follows from dominated convergence and the inequality < 1 — e Sy / S < Sy/s. Similarly, for y < and 
6 > 0, 



(1-5 (<%)) 



— ) 



5\0 



To show (2.11 i), note that 



lim q (dy) 1 ^ = lim (1 — 5 



i/<5 



{y>o} 



,-|»|G n 



{l/<0}- 



□ 



Then, (2.11 ii) follows from ( |2 . 1 If i) by taking logarithms on both sides. 

Following a density transformation construction as given in Sato [32] (see Definition 33.4 and Example 33.14 there) 
and using the martingale condition Ee^' = 1, we define a probability measure P* on (fi, J 7 ) via 



dP* 



(t > 0); 



(2.12) 



i.e., P* (B) — E (e Xt 1b) , for any Be J f and t > 0. The measure P* can be interpreted as the martingale measure when 
using the stock price as the numeraire. Under this probability measure, (X t )t>o has the representation (erW t * +L^)t>o 
where, under P*, W* := {W^) t >o is a Wiener process and L* := (L% ) t >o is a Levy process with triplet (0, b* , v*) given 
by 

v*{dx) := e x v(dx) = e x s(x)dx, b* := b + x (e x - 1) s(x)dx + a 2 , 

J\x\<l 

which is moreover independent of L* . Hereafter, we set 

<f (a) = q{x)e x , s*{x) = \x\- Y ~ l q* (x), q*(x) := q(x)e x . 



Remark 2.3. As it turns out, under ( 2.4 1 (equivalently, under (2.5)), v* corresponds to a tempered stable Levy 
measure. Indeed, setting T(s) := (- — l) -1 and, for x > 0, 



g*(x) := e x g(a;) = e x / e -5s r p(ds) = / e'^- 1 ) s Y p{ds) = / e --uV(du) 



with p* (cZu) | (q,oo) : ~ l(0,oo ) (u)(l + u)- y /0oT -1 (du), while, for a; < 0, 



/0 /-O /■() 

e~ x/s \s\ Y p(ds) = e-<i- 1 ){-s) Y p(ds)= e~ ^ {-u) Y p* {du) , 
CO </ — OO </ — OO 



(2.13) 



(2.14) 



with p*(du)\ ( _ oofi) := l ( _ 1)0 )(u)(l + «) y poT ^du). 



The following simple consequence of Lemma 2.2 will also be needed in the sequel: 
Lemma 2.4. Let 

M* := M - 1 = -q'(0 + ) - 1 and G* := G + 1 := <f(0") + 1, 



and Zei i/ie conditions (2.4), \2.6-i), and (2.6-ii) be satisfied. Then, M* > and, for any y € Mo; 



1/5 _ - S M*. 



(i) I™ 9* (<fy) ' = ' 



{y>o} + e 



-l^ G *l { , <0} , (ii) \imJ-\nq*(Sy) = -yM*l {y>a} -\y\G*l {y<0} . (2.15) 
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Proof. The first assertion is clear since, as it was shown above, (2.4) implies (2.5) and, therefore, using (2.9), 

p(ds) > 0. 



jr-*-i-_£V(i-i)**)-j[V(i-i 



Next, from Lemma 2.2 



1 



g (1 - r (<$»)) = I (1- ?(«tf)e ff ") = ^(1- + ^ (1- ^) ^ V(M - l)l y>0 + \y\(G + 1)1„ <0 . 



Using the above and proceeding as in the proof of Lemma 2.2 lead to (2.15). □ 

An important tool thereafter is to change probability measures from P* to another probability measure, say P, 
under which [L*)t>o is a stable Levy process and [W*)t>o is a Wiener process independent of L*. Concretely, let 



i>{dx):=\x\ Y 1 (C+l {x>Q } + C_l {x<0 }) dx, b = b* 



x(p — v*){dx). 



\x\<l 



Note that v is the Levy measure of a F-stable Levy process and that 

v{dx) = e (p ^v*(dx), with ip(x) := —hiq*(x) = — lnq(x) —x. 

Hence, by virtue of Theorem 33.1 in [32], there exists a probability measure P locally equivalent]^] to P* such that 
(X t ) t >Q is a Levy process with Levy triplet (<r 2 , b, v) under P. Throughout, E denotes expectation under P. 

Remark 2.5. By virtue of |421 Theorem 33.1], the previous change of probability measure is valid since, under our 
assumptions, the following condition is satisfied: 



(e^ )/2 -l)V(dx)= j (>(*V 2 -l) 2 ' e -^u{dx) 



< 00. 



Indeed, the integral above can be expressed as 

2 



a 







(e5^"*W_i) Ixl^dx + C. (ei^^-l) M^ 1 



dx, 



(2.16) 



as x \ 0. 



and, in light of (2.15), the integrand in the first integral is such that 

(e^™ -I)' ~i(lnr(*)) a ~ { ^x\ 

This shows that the first integral is finite on any interval (0,e). Outside any neighborhood of the origin, this integral 
is finite in view of (2.4). The second integral in (2.16) can be handled using again (2.15) and the fact that lug* (a;) < 
when x < 0. 

Denoting 7 := EXi = EL\, we recall that the centered process (Z t )t>o, defined by 

Z t := L* t - ty, (2.17) 

is a strictly ^-stable under P and, thus, is also self-similar; i.e., 

(t x l y ' Z ut ) u > = (Z u ) u >o, (2-18) 

for any t > 0. Let us now recall some standard results on stable Levy processes. Throughout, let pz denote the 
marginal density function of Z^ under P. It is well known (see, e.g., (14.37) in [15] and references therein) that 



Pz{v) ~ C±\v\ 



-Y-l 



as v — > ±00, respectively. 



(2.19) 



1 Equivalently, there exists a process (Ct)t>o such that P*(_B) = W(e Ut \g), for t > and B S Tt- 
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In particular, for any v > 0, as t — > 0, 

P(Z t > «) = P (Zj > t~ 1/y ~v) - «/([«, oo)) = t -y-v~ Y , (2.20) 

P(Z t < - v) = P < - t~ 1/Y vj ~ tP((-oo, -«]) = t^\v\- Y . (2.21) 
The following tail estimate can also be deduced from (2.19) (see Appendix [C| for its proof): 

P(|3t| >^ 1/y w) =P(|Z t | >«) <Kt|«r y , (2.22) 
for all v > and < £ < 1 and some absolute constant n < oo. 

Of use in what follows is the following representation of the density process (see Theorem 33.2 in [42] ) 



d¥* 



(2.23) 



with 



s<t 

AX 3 |> 



The process (Ut)t>o can be expressed in terms of the jump-measure N(dt, dx) :— #{(s, AX S ) G dtx dx} of the process 
(Xt)t>o and its compensated measure N(dt,dx) :— N(dt,dx) — v(dx)dt (under P). Indeed, 



U t = lim 



(ff tp(x)N(ds, dx) ~ t f (e^ x) - l) e-'^ x) v{dx)\ 

\Jq J\x\>s J\x\>e V ' ) 

t 



— lim 

£->0 



(p(x)N(ds, dx) + t / fe-vW - 1 + ip(x)) v{dx) 

IQ J\x\>e J\x\>e * ' / 

Finally, from the dehnition of Poisson integrals under compensated Poisson measures as in |23j (see Theorem 10.15 
therein), we get 

U t = U t + rit, 

with 

rt 



(2.24) 



U t := / tp(x)N(ds,dx), n:= / (e - ^ - 1 + <p(x)} u{dx) 



(2.25) 



Remark 2.6. Both Ut and 77 above are indeed well defined since from (2.4), (2.6-iii), (2.15 ii), and since ip(x) > for 
all x G M_ and Y G (1,2), 



kl<i 



- 1 + ^(x) |a;|- y - 1 da; < Ci / <^ 2 (x)|a ; r y - 1 da; < C 2 / |a;| 1 - y da; < 00, 



kl<i 



a<-i 



-<p(x) 



- l + <p(x) \x\- r ~ 1 dx < d / |<^(a;)||a;|- y - :L da; < C 2 / |a:|- r dx < 00, 



s<-i 



X<-1 



X>1 



e-vM -l + <p{x) \x\^ Y ^dx < C x / e lnq(x)+x x- Y - 1 dx+ / |-1 + <p(x)\ x^^dx < 00, 



X>1 



X>1 



for some constants C\,C% < 00. 



Finally, let us also note the decomposition of the process X in terms of the compensated measure N(dt, dx) 
N(dt,dx) - v{dx)dt: 



X t = aW t * + Zt+tj^ aW* 



xN(ds, dx) + 1"/. 



(2.26) 



The representation (2.26) can be deduced from the Levy-Ito decomposition of the process (X t )t>o (cf. Theorem 13.4 
and Corollary 13.7 in [23J or Theorem 19.2 in [42]). the stated Levy triplet (0, b, v) of X under P, and since b must be 
such that E(X t ) = jt, by construction. 
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3 The pure-jump tempered stable model 



In this section, we find the second-order asymptotic behavior for the at-the-money call option prices (1.2) in the pure- 
jump tempered stable process (i.e. a = and X t = LJ). The proofs of all the results of this section are deferred to 
Appendix [Aj Before stating our first result, we need to rewrite the call option price (1.2) in a suitable form. 

Lemma 3.1. With the probability measure P defined in (2.23) and the parameter 7 := E_Xj = EL|, 



t-±-^E(S t -So)+ = e -^ t 



E e 



-u t - 



_^ t l-l/Y 

The next two lemmas are crucial to obtain the main result of the section: 
Lemma 3.2. In the setting and under the assumptions of Section^ for any £ > 0, 

- i/Y {u t+Zt )\ = n e 



{t-l/YZ t >v} 



dv. 



(3.1) 



(i) limE(e- ?r 



ii) limE(e-«* 1/YUt 



t-i-0 



(3.2) 



where fj := T(-Y) (C+ M Y + C_ G Y ) and 77* = T(-Y) (C+{M*) Y + C_(G*) y ). In particular, both t~^ y ' (U t + Z t ) 
and t~ 1 / Y U t converge to a Y -stable distribution as t — > 0. 

Lemma 3.3. In the setting and under the assumptions of Section^ the following two assertions hold true: 



1. For any v > 0, 



lim 

t->o 



U t >v)= t-{x--\nq(x)>v}v{dx) 



(3.3) 



*-{ — \n q(x)>v}% 



-Y-l 



dx + C- 



L{ — x— \n q(x)>v} \X\ 



-Y-l 



dx. 



2. There exist constants k < 00 and t > such that 

1; 



t 

for any < t < to and v > 0. 



(i) -P [Ut > v < kv~ Y , (ii) -P [Zl + Ut > v < kv 



(3.4) 



We are now in a position to establish the main result of this section. The following result gives the second-order 
asymptotic behavior of at-the-money call option prices under an exponential tempered stable Levy model. 

Theorem 3.4. Under the exponential tempered stable model without Brownian component (a = 0), 



lim t Y 1 
t->o 



t-r E(S t -S ) + -E{Z+) ) = i? + t ? + 7 P(Z 1 > 0) , 
<~>o 



where n is given as in (2.25), 7 := EL\, and 

d = C+ / (1 - q{x) + \nq{x)) x~ y - x dx + C_ / (1 - e x q(x) + bxq(x) + x) \x\ 



-Y-l 



dx. 



(3.5) 



(3.6) 



Remark 3.5. The leading term E(Z ] f ) can be explicitly computed via the absolute first moment of Zy. 

E\Z X \ = 2E(Z+) - E{Z X ) = 2E(Z+). 

Using (1.2.13) in [41], it turns out that 



E{zt) = \nzi\ = -[c + + c.y rli- 



Y 



I 1 / C+ - c_ 

x cos — arctan ^ ^— tan 



C+-C- 
C+ + C- 

Yn 



tan 




(3.7) 



where C± = C±T(— Y)\ cos(iY"7r)|. In the symmetric case (i.e. C_ = C+ = C), (3.7) takes the form (1.6). 
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Remark 3.6. From (3.5), it follows that the short-time second-order asymptotic behavior of the ATM call option 



price (1.2 1 has the general form: 



1 



E(5 t -5 ) + =dit* +d 2 t + o(t), (*->0), 



(3.8) 



with d\ = E(Z+ ) given as in (3.7) and d 2 = ■& + i] + jF(Zx > 0). In words, the first-order term synthesizes only 



the information on the tail-heaviness index Y and the overall jump-intensity parameter, while the second-order term 
incorporates also the information on the relative intensities of negative and positive jumps as measured by C+ and C_. 



We now proceed to study the asymptotic behavior of the corresponding Black-Scholes implied volatilities. Through- 
out, let a(t) denote the ATM Black-Scholes implied volatility at maturity t with zero interest rates and dividend yield. 
The following result gives the asymptotic behavior of a(t) as t — > 0. 



Proposition 3.7. Under the exponential tempered stable model (1.3) without Brownian component, the implied volatil- 
ity a has the following small-time behavior: 



a(t)=aitY i + <7 2 t 2 + o(t*), 



t->0, 



(3.9) 



whe 



ax := V27rE(Z+), 

a 2 := V2^ (tf + r] + 7P(Zi > 0) 



(3.10) 
(3.11) 



and where & and rj are given by (3.6) and (2.25), respectively 



4 Tempered stable model with Brownian component 



In this part, we consider an exponential tempered stable model with non-zero Brownian component. Concretely, 
throughout, (X t )t>o is a Levy process with triplet (a 2 ,b, v) as introduced in Section [2] and u ^ 0. In that case, it is 
well-know that 

limE* I exp | iu—^ ) ) = exp 

*-»«> V V VtJJ 



2 2 

-a u 



and, therefore, X t /Vt converges weakly to the centered Gaussian distribution with variance a 2 . Equivalently, recalling 
that under P*, (Wj*)t>o is a standard Brownian motion, it follows that 



lim I 



ft 



> X = 



•{aW^ > x). 



(4.1) 



The first order asymptotic behavior for the ATM European call options in this mixed model is obtained in [33] using 
Fourier methods. We give next the second-order correction term for the at-the-money European call option price. As 
before, we change the prob ability measure P* to P so that X t — tj + aW? + Z t , with (Z t )t>o a strictly y-stable Levy 
process under P (see (2.17)). Recall also that, under both P* and P, W* is still a standard Brownian motion. We will 



also make use of the decompositions ( 2.24 )-( 2.26). The proof of the result below is presented in Appendix IBl 



Theorem 4.1. In the setting and under the conditions of Section^ the ATM European call option price is such that: 

C-+C+ 



t-S-0 



1 

So 



'(/ : 4- nS t - S )+ - aE*(W?) + j = 2Y ^_ + a 1 '* E* flW x * 



1-Y\ 



(4.2) 



Remark 4.2. The (1 — F)-centered moment of a standard normal distribution is given by (see, e.g., (25.6) in [42]): 

E* (|IF 1 *| 1 - r ) 



2!- y / Y 

1- — 

\/tt V 2 
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Thus, the second-order asymptotic behavior of the ATM call option price (1.2) in short-time takes the form 



■ E(S t - S )+ = dit* + d 2 t^ + o (t^) , (t -> 0), 



(4.3) 



with 



,i-y 



di 



do = 



r l- 



(4.4) 



Y\ (C_ + C 4 , 
/2tt ' ■ V 2 J 2F(y - 1) 

Intuitively, the first-order term only synthesizes the information about the continuous volatility parameter a, while 
the second-order term incorporates also the information on the tail index parameter Y and the overall jump-intensity, 
which is measured by the parameter C := C_ + C+. However, these two-terms do not reflect the individual intensities 
of negative or positive jumps as measured by the values of C_ and C+. This fact suggests that it could be necessary 
to develop a third-order approximation. For more information on the latter in the CGMY model, we refer the reader 

to [ig. 

The next proposition gives the small-time asymptotic behavior for the ATM Black-Scholes implied volatility under 
the exponential tempered stable model with non-zero Brownian component. Unlike the pure-jump case, we can only 

In fact, the first order term of the ATM call option price under 



derive the first order asymptotics using Theorem 4.1 



the generalized tempered stable model is the same as the one under the Black-Scholes model. The second order term 
requires higher order asymptotics of the ATM call option price. The proof is deferred to Appendix |B) 



Proposition 4.3. Under the exponential tempered stable model (1.3) with non-zero Brownian component, the implied 
volatility a has the following small-time behavior: 



a(t) 



Y(Y-1) 



-r i 



t L 



oU 1 ' 



t -> 0. 



(4.5) 



Remark 4.4. A careful analysis of the proofs of the main Theorems |3 . 4| and |4. 1 1 reveals that the results therein remain 



true provided that the assertions of Lemmas 3.2 and 3.3 hold true and that the two changes of proba bility meas ures 
of Section [2] first from P into P* and then from P* into P^ are valiq^] The first change is valid if |2j| and 



are satisfied, while, as seen in Remark 2.5 the existence of P is guaranteed if (2.10) (which in turn implies (2.15)) is 



satisfied with M > 1 and \nq*(x) < 0, for all x < 0. In particular, these two conditions are satified if ( |1.11[ ) and 



q(x)<C+, q(-x)<C-, 



(x > 0), 



(4.6) 



(2.6 in) and that 



hold true. On the other hand, the proof of Lemma 3.2 will remain valid if, in addition to (1.11) and (4.6), we assume 

q(x) < C+e- x , (x > 0). 



Finally, Lemma 3.3 does not require any additional assumption. 



5 Asymptotics for ATM option prices under a CGMY model 



In this section, we specialize the second-order asymptotic expansions of the two previous sections to the CGMY model. 
The result presented here were first reported in [16] . Recall that under the CGMY model, the Levy measure of the 
pure-jump component (L t )t>o is given by 



v(dx) = \x\ q(x)dx = \x\ 



-Y-l 



(Ce 



-Mi • 



l {x>0} 



Ce 



Gx 



L {a<0} 



dx, 



(5.1) 



with corresponding parameters C, G, M > and Y £ (1,2). Then, the characteristic function of the log-return process 
X t = aW t + L t takes the form 



ip t (u) =E(e ruXt ) =cxp (t 



ICU 



+ CT(-Y) ((M - iu) Y + (G + iuf 



M — G 



(5.2) 



2 Indeed, the first transformation enables us to use the representation of Lemma |3.l| while the second one enables us to exploit the 
self-similarity of the process (Zt)t>o an( i its more explicit distributional behavior in short-time. 
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for a constant c£l. The martingale condition (2.1 1 implies that M > 1 and 



-CT(-Y) ({M - if + (G + 1) Y - M Y - G Y ) 



2 ' 



(5.3) 



see, e.g., Proposition 4.2 in 03] . In particular, the center 7 := EXi = ELi of X and the parameter b of X (relative to 
the truncation function irl.r| x |<n) are given by 

7 = c- CYY{-Y){M Y ~ l - G Y - 1 ), 



b = c 



xv{dx) - CYT(-Y)(M 



, r Y-l Q Y ~ X 



(5.4) 
(5.5) 



|x|>l 



Under the share measure P* introduced in Section[2j (X t )t>o has Levy triplet (b* , (a*) 2 ,i>*) given by 



a* := a, b* := c + er 2 



xv*{dx) — CYT{—Y){{M*) — (G*) ), u*{dx) := e x u{dx), (5.6) 



|x|>l 



with M* — M — l and G* = G + l. Finally, under the probability measure P, the centered process (Z t ) t >o is symmetric, 
and its center 7 := EXi = EL\ is given by (see [16] for the detailed computation) 



7 = 6 + 



{|x|>l} 



xv{dx) = b* + I x(y — v*)(dx) + 



|x|<l 



= -CT(-Y)((M -If +{G + l) Y -M Y -G 



{\x\>i} 

Y' 



xv{dx) 



(5.7) 



We also need the value of 77 defined in (2.25), which, under the CGMY model, is now given by 



r) = C (e 
Jo+ 



-M*x 



l + M*x)x- y - i dx + C (e° X -1-G*x)\x\ 



-Y-l 



dx 



= CT{-Y)({M*) Y + (G*) y ), (5.8) 

where, for the last equality above, we used the analytic continuation of (14.19) in [32] , 

We are ready to explicitly write the second-order e xpansi ons. First, let us compute the term d of (3.6). To this 
end, it is convenient to use the representation given in (A. 15 1 below noting that q(x) := e~ Mx l{ x> oj + e Ux l{ x< o} in 
the CGMY case: 

pOO f 

-Y-l, 



d = C I (e v - 1) / l{ x --\nq(x)>v}\x\ dxdv 



CI (e v -1) I l {Mx > v} x Y - 1 dxdv + C I (e v - 1) I l { _ x _ Gx > v} (-x) Y 'dxdv 



\-Y-l, 



c 

Y 



:(M y + (G + l) y ) j {e- v -l)v- Y dv 

= -CT(-Y) ((M* + 1) Y + (G*) Y ) , 
where in the last equality we used the following identity (see p. 84 in |4"2"]): 

/>oo 

/ (e~ v — l)v~ Y dv = T(l — Y) = —YT(—Y). 
Jo 

Using ( 5.7 1-( 5.9 ), it follows that in the pure-jump CGMY model, 

7 Gr(-Y) 



(5.9) 



d 2 := 1? + rj + ± 



((M - i f - M Y - (G + i f + G Y ) 



while in the general CGMY model with non-zero Brownian component, G+ = G_ = G, and thus 



dn : = 



Ca 1 - Y 2 1 - Y 
Y(Y-l)v^ J 



r 1- 



Y 
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6 Numerical examples 

In this part, we assess the performance of the previous approximations through a detailed numerical analysis for a 
CGMY model. 



6.1 The numerical methods 



Let us first select a suitable numerical method to compute the ATM option prices by considering two methods: Inverse 
Fourier Transform (IFT) and Monte Carlo (MC). 

Before introducing the IFT method, let us set some notations. The characteristic function corresponding to the 
Black-Scholes model with volatility £ is given by 



ip t (v) = cxp -- 



(, 2 



The corresponding call option price at the log-moneyness k — \og(So/K) under the Black-Scholes model with volatility 
£ is denoted by Cg S (k) ; that is, 



C% s (k) = S e- rt E(e 



rt V (Jr-T, 2 /2)t+ZW t _ 



Let use also recall that the characteristic function under the generalized CGMY model is denoted by <p t (see (5.2)) 



and let us denote the corresponding call option price at log-moneyness k by C(k). The IFT method is based on the 
following inversion formula (see Section 11.1.3 in [8]): 

1 



z T (k) :=C(k)-C^ s (k) = 



2tt 



-ivk 



( T (v)dv, 



(6.1) 



where 

W (v — i) — w BS ' S (v — i) 
iv(i + IV) 

In our case, we fix r = and, since we are only interested in ATM option prices, we set k = 0. In order to compute 



(6.2) 



numerically the integral in (6.1), we use Simpson's rule 

z T (0) := 



1 

27 



C T (v)dv = A 



p-i 

771=0 



M 



C T (v 



Mh 



with A = Q/(P - 1), v* M = -Q/2 + mA, and w Q = 1/2, w 2 e-i = 4/3, and w 2 t = 2/3, for £ = 1, . . . , P/2. 

We also consider a Monte Carlo method based on the risk-neutr al opt ion price representation under the p robab ility 



measure P. Under this probability measure and using the notation (2.25) as well as the relations (2.24) and (2.26), we 
have: 



E(e 



1) 



E" 



x 7 



1 



E 



T 



= E ( e ~ M 'U}+G*U-- V T f 1 _ e -V+-U--T^-aW T 



the variables and 



U T are 



which can be easily computed by Monte Carlo method using the fact that, undcr 
independent Y-stable random variables with scale, skewness, and location parameters TC\ cos(ttY/2)\T(—Y), 1, and 
respectively. Standard simulation methods are available to generate stable random variables. We consider the following 
set of parameters 

C = 0.5, G = 2, M = 3.6, Y = 1.5. 

Figure [T] compares the first- and second-order approximations as given in Remarks |3.6| and |4.2| to the prices based 
on the Inverse Fourier Transform (IFT-based price) and the Monte Carlo method (MC-based price) under both the 
pure-jump case and the generalized CGMY case with a = 0.4. For the MC-based price, we use 100, 000 simulations, 
while for the IFT-based met hod , we use P = 2 14 and Q = 800. As it can be seen, it is not easy to integrate numerically 



the characteristic function (6.2) since T is quite small and, therefore, the characteristic functions ipx and (p T 
quite flat. The Monte Carlo method turns out to be much more accurate and faster. 
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Pure-Jump CGMY Model (C=0.5,G=2,M=3.6,Y=1.5) 



General CGMY Model (o=0.4,C=0.5,G=2,M=3.6,Y=1.5 



j35 



□ 035 




IFT-based Method 



IFT-based Method 

MC-based Method 

1st-order Approx. 

2nd-order Approx. 




0.5 



1 1.5 2 2.5 3 



3 5 



0.5 



1 1.5 2 2.5 3 

Time-to-maturity, T (in years) 



3.5 



Time-to-maturity, T (in years) 



x 10 



Figure 1: Comparisons of ATM call option prices for two methods (Inverse Fourier Transform and Monte-Carlo 
method) with the first- and second-order approximations. MC-based price is based on 100, 000 simulations while the 
IFT-based method is based on the parameter values P = 2 14 and Q = 800. The parameter a in the generalized CGMY 
model is set to be 0.1. 

6.2 Results for different parameter settings 

Here, we investigate the performance of the approximations for different settings of parameters: 

1. Figure [2] compares the 1st- and 2nd-order approximations with the MC prices for different values of C, fixing the 
values of all the other parameters. In the pure-jump case, the 2nd order approximation is significantly better for 
moderately small values of C, but for larger values of C, this is not the case unless T is extremely small. For a 
nonzero continuous component, the 1st order approximation is extremely bad as it only takes into account the 
parameter a. 

2. Figure [3] compares the 1st- and 2nd-order approximations with the MC prices for different values of Y, fixing 
the values of all the other parameters. In both cases, the 2nd order approximation is significantly better for 
values of Y around 1.5, which is consistent with the observation that \d%\ — > oo as Y — >• 1 or Y — s- 2. For a 
nonzero continuous component, the 1st order approximation is again extremely bad as compared to the 2nd order 
approximation . 

3. In the left panel of Figure [4j we analyze the effect of the relative intensities of negative jumps compared to 
positive jumps in the pure-jump CGMY case. That is, we fix the values M to be 4 and consider different values 
for G. As expected, since the first order approximation does not take into account this information, the 2nd-ordcr 
approximation performs significantly better. 

4. In the right panel of Figure|4j we analyze the effect of the volatility of the continuous component in the generalized 
CGMY case. The 2nd order approximation is, in general, much better than the lst-order approximation and, 
interestingly enough, the quality of the 2nd order approximations improves as the values of a increases. In fact, 
it seems that the 2nd-order approximation and the MC prices collapse to a steady curve as a increases. 
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Figure 2: Comparisons of ATM call option prices with the short-time approximations for different values of the jump 
intensity parameter C . 




Figure 3: Comparisons of ATM call option prices with the short-time approximations for different values of the 
tail-heaviness parameter Y . 
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Pure-Jump CGMY (M=4, Y=1.5, C=0.5) 



General CGMY (C=1, G=2, M=2, Y=1.5) 



— ^MC Price (G=M) 

- © - 1st Approx. (all cases) 

-^0- 2nd Approx. (G=M) 
— B— MC Price (G=M/4) 

-□- 2nd Approx. (G=M/4) 
MC Price (G=4*M) 

-z- 2nd Approx. (G=4*M) 




— e— MC Price (ct=0.1 ) 

- e - 1st Approx. (o=0. 1) 

- e - 2nd Approx. (o=0.1 ) 
— * — MC Price (tr=0.3) 

- » - 1st Approx. (o=0.3) 

- * - 2nd Approx. (cr=0.3) 

MC Price (o=0.5) 
1st Approx. (o=0.5) 
2nd Approx. (o=0.5) 
—i— MC Price (o=0.8) 

- * ~ 1st Approx. (o=0.8) 
* - 2nd Approx. (o=0.8) 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

Time-to-maturity, T (in years) 



_--*--"'*""" „ 1=1 -1-1- -n- B — Q -O- Q ~ 
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Figure 4: Comparisons of ATM call option prices with the short-time approximations for different values of G or 
and different values of the volatility parameter a\ 

A Proofs of Section [3^ Pure-jump tempered stable model 

For simplicity, unless otherwise stated, throughout this section, we fix So = 1. 
Proof of Lemma 13. 11 



From (1.9), we have 



t-^ Y E(S t - S ) + - t-V Y F*(X t >E)= t~ x / y / 

Jo 



e~ x F*(X t > x)dx 



: (t~ 1/Y X t > u)du. 



Next, using the change of probability measure ( |2.23[ ), 

v*{r 1 ' Y x t >u) = e* (i {t -vY Xt >u } ) = e (e- u n {t - 1/YXt > u} ) , 



and, moreover, from (2.171 and (2.24), 



\t- x ' Y X t >u) = e-"*E (e- Ut \ { 



t~ 1 / Y Z t >u~yt 1 - 1 / Y } 



Then, 



t- 1 / y E(5 i - S )+ = J e- tl,Yu -^ E (tr u < l {t - 

and, changing variables (y = u — 7i 1_1 / y ), the result follows. 
Proof of Lemma 13.21 



1 /yz t > u -7ti-v>'} ) du, 



(A 



(A 
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Fix Q t '■— Ut + Z t . From the usual moment generating function formula for Poisson integrals, the change of variables 
y = t~ x l^ x, and recalling that ip(x) + x = — \nq(x), 



E 



(e-e*- 1/y O*)=E( 



= E exp -£t 



exp 
exp 



{f(x) + x) N{ds, dx) 
t j (e£*- 1/Y to-fa) - l -fr^ln^x)) v(dx) 
C + ^ (e^ 1>Y ^ 1/Y y) l-^^q^y)) 



y- Y - x dy 



+ C-f (f- 1,Y M* ,Y *) l-^ Y \nq{t^ Y y)) ly^dy} . (A.3) 

Let S(u) := e" — 1 — u and _B := sup | lng(u>)|/|u>|, which is necessarily finite in light of (2.6 -in) and Now, 
using that \nq(u) < 0, u € Rq, and that |H(u)| < X (u 2 A for w < and some constant K < oo, 



£i v In q ( t v y 



\y\ 



-Y-l 



< K 

< 



((f»B) 2 A|fyB|)| 1 ,|- r - 1 
x((£5) 2 V|6B|) (y 2 A|y|) 



Y-l 



which is integrable. Therefore, one can pass the limit inside the integrals in (|A.3|) and, using (|2.1l|), 
= exp < C_j 



limE (e~ & 
t->o \ 



-i/y. 



= exp{ C. 



|c+ J™ (e-*t M - 1 + y^M) y- Y - x dy + C_ J (e-l^ G - 1 + \y\£G) ly^dy 

(c + J™ (e- uM - 1 + uM) u- y - l du + C_ J° (V^G _ 1 + | U | G ) \u\- Y - l dv^j £ y J. 



Finally, using the analytic continuation of the representation (14.19) in [42], we show that the last expression is of the 
form exp(ry*£ y ), with rj* given as in the statement of the Lemma. 

For (|3.2|4i), we proceed as above to get 



E e 



y y 



.-Y-l 



dy 



+ C- f EUt-^ Y \nq*(t^ Y y))\y\- Y - x dy 

J —OO 

Now, for all y e Mo, 

er 1 ^ in,-* (iV Y v) = er 1/y in?- (t^ Y y ) +^ = ^(^1^ {^ Y y 

Indeed, this is clear when y < 0, since lng(u) < 0, while for y > 0, (A. 5 1 follows from 

hiq~(x) 



1 < 0. 



(A.4) 



(A.5) 



sup 

a;>0 X 



< -1, 



since, from (2.5), for all x > 0, 

e x q{x) = e x f e~ x ' ' s s Y p(ds) = ( e~ x ^- 1 h Y p(ds) < [ s Y ' p{ds) = 1. 
Jo Jo Jo 

Using (A.5), we can proceed as above to justify passing the limit into the integrals in ( A.4[ ). Next, using Lemma 2.4 



limE 

t->o 



= exp < C- 



,-vW _ 1 + y £M*) y^^dy + C 



'^ G ' -l + \y\^G*)\y\ 



\-Y-l 



exp i C. 



(c+ J™ (e~ uM * - 1 + uM*) u- Y - l du + C- J° (e-l"l G * - 1 + |u|G*) \u\~ 



Y-l 



dy 



du C 
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Finally, (3.2 ii) follows once more from the analytic continuation of (14.19) in [22]. □ 
Proof of Lemma 13.31 

We use a small/large jump type decomposition. Concretely, fix e > and let 

z[ e) := f f xN(ds,dx), Z[ e) := Z t - z[ e) . 

JO J\x\>e 

Denote respectively by (iVj )t>o and )i>i the counting process and sizes of the jumps of (Z^ )t>o> so that 
Zf = £t ■ Next, defined the corresponding processes for U: 

uf := [ I <p{x)N{ds, dx) = - y £ i mf M e) ) , := U t - uf. 

JO J\x\>s i=1 K 7 

Let A e := EN[ £) be the jumps intensity of We prove the validity of the two assertions in two steps: 
(1) By conditioning on iV t , 

r l V ( Z+ + U t >v) = t^P ( ( Zl £) ) + + > v ) e-^ 1 



+ e- x ^\ £ V ((Z t (£) + + + C/ t (£) - lnf (^ e) ) > v\ + 0(t). 



The first term above can be made 0{t) by taking < e < £o, for some £o > (see, e.g., [3H1 Lemma 3.2] or JJH Section 
26]). For the second term, first note that for some £q > 0, small enough, 



a (v) := A £ P - mg* (f^) > A 

1{x+-\nq*(x)>v}v{dx) = I 1{ X - _ ln q( x )>v} v(dx) 



poo />0 
= C+ / l{-l nij ( x )> v }X^ Y ^ 1 (dx) + C- / l{-x-lnq(x)>v}\x\~ Y ~ 1 dx, 
JO J -oc 

for all < e < £o. Also, since 

F v , e (z,u) :=P^(z + ^) + +u-\nq* (^ s) ) > J) 
is continuous at (z, u) = (0, 0), for any fixed < £ < £o, the function 

At(v) ■■= ^ ( Z t + U t >v^)~ a (v) 

is such that 

lim A t (v) = X E lim (^(z ( t e) + ^ e) ) + + - lnf (^ e) ) > vj - A e P - ^<f > v 

= \ e limE (F v , e (zf\ C/ t (e) ) - F Ol8 (0, 0)) 
= A £ E (limF v , e (zf.U^) - F„, e (0,0)) 

where we had used the dominated convergence theorem to obtain the last equality. 



19 



(2) Throughout, k > denotes a generic constant that may vary from line to line. First, note that (2.22) implies that 



-^(zt >v)< kv- Y , (A.6) 



for any < t < 1 and v > and some constants k < oo. So, it suffices to show the analog inequality for U. 
Let C :— C+ + C_. Using the decompositions in the previous step (1) with e — v/4, 



ip (u? > v/2) < -p (n^ * o) = \ (1 



e~ Aet ) < A e < C / |x|- r - 1 da: < 

|»|>o/4 



for some k < oo. Then, 



EU^^t lnq*(x)D(dx) =t lnq(x)v(dx) + tC- / x\x\~ Y - 1 dx + tC + I x^l'^dx, 

|a;| >e J\x\>e J x<-e J x>e 

and, since q[x) < 1, for all x, EZ7j < iC+ / 2 . >£ a;|a;| _y_1 (ia: = C+(y — \)~ 1 4 Y ~ 1 tv 1 ~ Y . Thus, whenever t and v are 
such that C+(F - l)- 1 4 Y - 1 tv 1 - Y < v/4, we have 

P(C/ t (£) > u/2) < V(ul E) - E L> t (e) + E C/ t (e) > u/2) < P(C/ t (e) - E £/ t (e) > v/4). 
Next, using a concentration inequality for centered random variables (see, e.g., |22j . Corollary 1), 



-2 



P(^ (£) > «/2) < P(# e) - > v/4) < e *H* + ^) los ( 1+ ^J < ( ^Yl) " t i < 

\ ev J ir 

where V/ := Var(J7 t ( - £ - 1 ) and in the last inequality we recall that e = v/4. Now, 



V 2 = Var(C/ t {£) ) = / {\nq*(x)) 2 D(dx) < C / x 2 \x\- Y - l dx < kv 2 - Y 

J{\x\<e] J{\x\<e} 



iii) and 



for some constants k < oo, where above, we set C — supu^g | lng*(x)|/|x|, which is finite in light of (2.6 
( |2.15f ii). Therefore, whenever C+(Y - l)- 1 4 Y - 1 tv 1 ~ Y < v/4 (or equivalently, C+(Y - l)-H Y tv~ Y < 1), 

\f(uP>v/2) < &V ^<kv- Y . 

Moreover, for any t > and v > 0, 

^P([7 t (e) > v/2) = ip(c7, (e) > «/2)l {c+(y _ 1) -i 4 y to -v< 1} + ip(L^ > «/2)l {c+(r _ 1) -i 4 v to -v >1} 
< + ^±(Y - l)-H Y tv- Y < k'v' Y , 

for some constant k . Combining the previous estimates, we finally have 

-V(u t > v) < -$(u* > v/2) + jV(u t ie) > v/2) < kv~ Y , (A.7) 
for all v > and t > and some constant k < oo. □ 
Proof of Theorem 13.41 



For simplicity, we first treat the case 7 = so that, in light of Lemma 3.1 



/■OO _ 

r'^EiSt - S ) + = e~* / e' tl/Yv E (e^l^-w^j) dv. (A.8) 
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The general case is resolved in Lemma |A.1| below. Let 

D{t) := t-V Y -E(S t - S a )+ - E(Z+), 

which can be written as 

D(t) := jf e- tl/Yv E (e-^l^^^^j dv - E{Z+) 
+ (e-"« - 1)E(Z+) 

+ (e- 7>t - 1) f jH e- 4l/ ^E (e-^ l^-vy^j) - E(Z+)\ (A.9) 
=:Di(t)+I>2(t)+Ds(t). 

We will show that 

t^^Dxit) — ><& + r], as t->0, (A.10) 
for a certain constant while it is clear that D^(t) = o(Di(t)) and that t 1 ' Y ~ 1 D2{t) = o(l), as t — > 0. First, 

D x (t) = E ( e- Ut I e~* v dv - E(Z+) 



t~ x l Y ri(e- & *)-E(e-( Di - 

i/y ^t_g( e -(^ + ) 



E(Z+) 

i" 1 ^ hp -E{e-V t+ *< ) ) ) - E(Z+). (A.ll) 

Next, by the self-similarity of (Z t ) t >o under P (see \2.18\ ), E(Z+) = t~ l / Y E{Z+) and, also using that E(U t ) = 0, 
, / e ^_i 1 - E (e-^ + +^)) - E (Z+ + U t 

= + t E I / - ^ dwl {^ + +^>o} ) - i E ( j z++ ~ t " !) dvl {zt + u t <o } ) 

= + - / (e- u - 1) P(Z+ + U t > v)dv - (e v - 1) P(Z+ + C7 t < (A.12) 

t t Jo t Jo 

=: D u (t) + D 12 (t) + D 13 (t). 



Clearly, 



Dn(i)->77, as t->0. (A.13) 



Let Qt := Z t + E/t and note that Qt < Zf + Ut- Then, for D\j,{t), using that e v — 1 < ye y , y > 0, Markov's inequality, 
and since 1 < Y < 2, 

/>00 

o < Di, 8 (t) < t 1 ^- 1 y (e fl/y " - 1) p (<- 1/r g t < -«) dti 

/•OO 

< t 2 / y - 1 J e tl/Yu uP (t-^ Y Q t < -uj du 

poo 

<t 2 / Y - 1 e( tI/r - 1 ^«d«E(e-*" 1/V «*) ^§ 0, (A.14) 



where in the last step we applied ( 3.2 i) . 
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Let us now deal with D12. First, from ( 3.4 ii) , one can apply dominated convergence and pass the limit inside the 
integrals so that 



(t) = r°(e-» - 1) lim (t-^(z+ +U t >vj) 



dv. 



From (3.3), wc then have 



UmD 12 (t)= / (e 11 - 1) / l{ x --ing(x)>'u}i > (^)^ := 1?. 

•/ «/ M.n 



(A.15) 



Combining (A. 12), (A. 13), (A. 14 1, and (A.15), it follows that 



\imt 1 ' Y - 1 D 1 {t)= (e~ v -l) l {x --inq(x)>v}Hdx)dv + ri. 

JO JRn 



Finally, the expression in (3.6) follows recalling that v(dx) = \x\ Y 1 (C+l{ a;> o} + C_l{ ;E< o}) dx and applying Fubini's 
theorem to the right-hand side of (A.15): 



C+ \ (e v -l) l{-\nq(x)>v}X dxdv + C- / (e v - 1) / l{- x -i n q(x)>v}\%\ dxdv 

JQ JO J0 J-00 

/>oo p— lnq(x) pO r—x—lnq(x) 

= C+ / (e~" - l)dvx~ Y - l dx + C_ / / (e"" - ^dw^r^- 1 ^ 

Jo Jo J -co JO 

= C+ f (l - e ln «^ + lng(.x)) sT^da; + C_ / (l - e^ 1 "^^ + x + \nq(x)^j {x^^dx. 



This concludes the proof. 



Lemma A.l. // 7 ^ m £/ien 



limt^- 1 t-^— E(5 t -S' )+-5 E(Z+) = tf + 7/ + 7 P(Z 1 >0). 
t— >o \ Oq 



□ 



(A.16) 



Proof. Without loss of generality, fix So = 1 and also assume that 7 > (the case 7 < being similar). Using (3.1 ), 
limi 1 /^ 1 (t-V Y K(S t -S ) + -E(Z+)) =lim< 1 /^-i f e -W+^ e -* 1/y *g (e"^ l {t -i/v Zt >„ } ) <fo - E(Z+) 

+ limtV^i e -(7 + ,) t f° e-^^re-^l^vr^j 
=:5 u (t)+5i 2 (t). 



As in the proof of (A. 10) 



limZ)ii(i) = $ + 77. 



(A.17) 



For Z?i2(i), changing variables to u = t 1 ' Y 1 and probability measure to P*, we have 

5, 2 ( t ,= e -./_; e -E( e -^.«» 1{z , 2M ),„ . «-*£ e -r(r- Z , > 

Next, recall that (L*)t>o and, therefore, also (Zt)t>o, are tempered stable Levy processes under P*, as shown by 
( 2.13 )-( 2.14 ). Then, since Y € (1,2), t~ x / Y Z t converges in distribution to a F-stable random variable Z under P* with 
center (or mean) and spectral measure 



tr(+l)= / u Y p{du) = s Y p(ds) = C+, er(-l) 
Jo Jo 



\u\ Y p(du) 



\s\ Y p(ds)=C-, 
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see Theorem 3.1 in [35] (see also Proposition 1 in [37]). Hence, the distribution of Z (under P*) is the same as the 
distribution of Z\ under P. Thus, Slutsky's lemma implies that t~ x l Y Z t — t 1 ^ 1 ^ u — Z and, in particular, 



lim P* 



(t- 1/Y Z t - t x - 1/Y u > o) = P* [z > o) = P (Z 1 > 0) 

Finally, by the dominated convergence theorem, 



limD 12 (t) = 7P(Zi > 0) . 



(A.18) 

□ 



The small-time asymptotic behavior of the ATM call option price CBs(t,a) at maturity t under the Black-Scholcs 
model with volatility a and zero interest rates is (recall that Sq = 1) 



Combining ( |A.18[ ) with ( |A.17| ) leads to ( |A.16[ ). 
Proof of Proposition |3.7[ 



C B s(t,a) 



/2tt 



24^2^ 



(A.19) 



see, e.g., [T5] Corollary 3.4]. To derive the small-time asymptotics for the implied volatility, we need a result analogous 
to (A.19) when a is replaced by ct(t). The following representation taken from [351 Lemma 3.1] will be useful, 



C BS (t,<T)=F(*Vt) with F(6):= &(-)dv = -= 

JO VZTT Jo 



exp 



dv, 



together with the Taylor expansion for F at 9 — (see [3j)l Lemma 5.1]); i.e. 

1 „ 1 



F{9) = 



+ O(0°), 0->O. 



/2tt 24^2^ 

Then, since a(t) — > as t — > (see, e.g., [13J Proposition 5]), we conclude that 

HQ +1/2 _ "'ft) 3 ,3/2 



CBs(*,*(t)) 



'2tT 



24V2 



= t3/a + of(^(i)tVa)' 

7T V 



as t -> 0. 



(A.20) 



Returning to the proof of Proposition 3.7 by equating (13. 81) and (A.20) and comparing the first order terms 



/2vr 



t->0, 



and, therefore, 



<r(t) - v / 2VE(Z 1 + )^~2 := fTi^"5, t^O. 



(A.21) 



Next, set cr(i) = cr(t) — u\t^ 2. By comparing the first and second order terms in (3.8) with the first term in (A.20) 
(noting that the second order term in (A.20) is o(t)), 

Hence, a(t) — > as t — > 0, and moreover 

cr(i) - v / 27^ + ?7 + 7P(Zi > 0)) t 1/2 , (A.22) 

□ 



Combining (A.21) and (A.22) finishes the proof. 
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B Proofs of Section [4j Tempered stable with Brownian component 



Proof of Theorem 14. XI 



For simplicity, fix So = 1. Recalling that X t = uW^ + L* t under P*, and using (1.9), the self-similarity of W* , and the 
change of variable u = i _1 ' 2 x, 



R t := r 1/a E(S t - S )+ - aE*(W?)+ = J e^T* (aW? >u- t~^L*J du - J 



'{<rWl > u)du. 



Next, changing the probability measure to P, using that L* t = Z t + jt, Ut = Ut + Tjt, and the change of variable 
y = u — t 1 / 2 ') in the first integral above, lead to 



Rt = 



= e -(>7+7)< 



+ e -(iJ+7)t 



Vty-jt jjjj g-^t-^t] 



{ffW 1 *>y-t _ 9^t} 



rfy- / Ee 



{<tW*>u} 



du 



Ee~ Ut l, 



(B.l) 



,-jt-Vty 



l P*W>y)4 



Above, the last term is clearly 0(i 1//2 ) as £ — > 0, while the sec ond te rm can be shown to be asymptotically equiv alent 
to a term that is 0(£ 1//2 ) by arguments analogous to those of (A. 18). Thus, we only need to study the term in (B.l) 
that we denote A t and which can be written as: 



At = e-^E e~ Ut 



-Viy 



( 1 {rrW*: 



.> v -t-iz t }- 1 i* w ?*vy) dv )' 

where we had set r) :— r\ + 7. To study the asymptotic behavior of At, decompose it into the following three parts 

r<rW*+t~^z t 



A t = e'^E e~ Ut l 



{W* >0,aW* +t~ 2 Z t >0 



o-Viy 



dy 



{o<ffW«<-r2Z t } 



<jW* 



-Viy. 



+ e-'''E e~ Ut l, 



e "-"dy 

aW? +t~ h Z t 



Viy d , 



y 



{0<-aW*<t 2 Zt } J Q 

=: hit) - h{i) + I 3 (t). 
We analyze each of there three terms in the following three steps: 

Step 1. We first study the behavior of h(t). Since {Z t )t>o and (W*)t>o are independent, 

/ e -u t _ e -(U t +Zt) \ 



(B.2) 



h(t) = e- m E 1 



+t~2Z t >0} 



poo 

J%(t,y)e 



-fjt 



{Z t >-t2y} 

V_ 

-Viy e 2 ° 



V2ira 2 



dy 



(B.3) 
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Using the self-similarity of (Z t )t>o & n d since EZ t = 0, Ji(t,y) is then decomposed as: 
J 1 (t,y)=E\l 



e -u t _ e -(u t +z t ) 



r^Zt ) I +^~5e((-Zi)i i i , 

1 ' V U {-Z!>t2-Yy} 



~{z t >-tiy}\ 
=■ Jn(t,y) + Ji2(t,y). 

Let us first consider J\2{t,y). From (2.19l- |2~21 1, there exists a constant A > such that 

for any < t < 1 and y > (see Appendix [C] for the verification of this claim). Moreover, for any fixed y > 0, 

j* oo poo 

t^'^ 1 Ji2(t,y)=t^ + ^^^ I upz(-u)du = i¥~"^~5 / wpz{-ti~^w)dw. 

Jt2~Yy Jy 



Using (2.191, there exists < to < 1 such that 

t^-v -5 wpz (-ti-^ w ) < 2(C+ V C_)u;~ r , 



for any < t < to and w > y. Therefore, by the dominated convergence theorem, and in light of (2.19), we get 

;dy 



limt^-V^ / J 12 (t ) y)e- Viy — 



V2~t 



oo / />oo 



(lim£2 2p z (— 



\Jy 

oo / />oo 



\t-s-0 



7T(J 



For J n (t,y), 



= C -lo U '' 



Jn(t, y) = HE ( l { , t > 0} £ +Zt (e- - - Hi ( l { _ t ^ t < 0} " X ) 



~ Y dw 



e 2* 2 , C 
dy = 



V2 



Y - 1 



x-y 



V2 



fdy- 



dx 



= H /" (e~ x - l)Ti(t,a;, y)dx-t~ $ f (e~ x - ljT 2 (t,x,y)dx, 

where, for t > and y > 0, we set 

T t {t,x,y) := V(Z t >0,Ut<x<U t + Z t ), T 2 (t,x,y) :=¥(- tiy < Z t < 0,U t + Z t < x < U t ). 



By (3.4-ii), there exists k > 0, such that for any x > and < t < 1, 

Ix(t,a:,j/) < P(x < J7 t + Z t ) < P(x < £7 t + Z+) < ktx~ Y . 



Hence, 



OO /-OO 



< e-Vf* 



< e~ m kr 



£-1 



/0 JO 

/•oo />oo 



T 1 (t,x,y)dx 



e -Viy e -^ 
V2na 2 



dy 



(1 - e- x ) X - Y dx-;^=dy — > 0, 

/0 J0 V27TCT 2 
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as t — > 0, since Y > 1. Similarly, using ( 3.4 i) , there exists a constant k < oo such that 

T 2 (t,x,y) < ¥(U t >x)< ktx~ Y ', 



for any a; > and < i < 1, and, thus, as in (B.9), 



lim 

t-s-0 t 1 - 



— f\t POO POO /-i — x 



2 Jo JO 



(1 - e~ x ) 

Vt 



T 2 (t,x,y)dx 



e -s/ty e -f 



V2t 



-dy = 0. 



For x < 0, using (3.2-ii) and Markov's inequality, there exist < to < 1 an d k < oo such that 



Ti(t,x,y) < ¥(r 1/Y U t < t~rx) < E(e" 
for any < i < t . Therefore, 

0<t^ I (e~ x - l)T 1 (t,x,y)dx < kt~r 

J —oo 

Hence, by the dominated convergence theorem, 

e -» _ i p-Vty 



')e f 



ty 



1-tr 



< e-^i 



4-1 



-T 1 (t,x,y)dx- 



V2ir<r 2 



dy < kt 2 



1-tT J V2 



- y dy — > 0, 



as t -> 0, since 2/Y > 1 > (3 - F)/2, for 1 < Y~ < 2. Similarly, using plfi), for x < 0, 

r 2 (t,x,y)<P(cJi + Z t <a?) <E(e- t " / "(^ +z 0)e t " / " a: <Ke*" /ra , 



for any < t < t and some constant k. Therefore, as in (B.13I, 



lime-"*^- 1 f f ^ B n ^ T 2 \l..r. ,,),!. v 

Jo J—oo 



t-¥0 



Vt 



dy = 0. 



Combining (B.6), (B.9), (B.ll), (B.13) and (B.14|, we finally obtain 



Y C 

t->o 1W Y-l Jo 



e 2a 2 
V / 2ttct 2 



dy. 



Step 2. Next, we study the asymptotic behavior of I 2 (t). Using the independence of (Z t ) t > and (Wj*) 

1 - e-VtaW? \ 



I 2 (t)=e- f > t E e~ Ut l 



e""* / E 
Jo 

POO 

e~ m / E 
Jo 



{0<aW*<-t-2 Z t } ^fl 

1 — e~^ y 



(e~ Ut l r i Y 

V {z t <-t2 a >y 



e 20 



-11. 



V< %/27rcr 2 
1 - e"^" . 



- 1 \ 1 - p" ^» 

• e -*y P(^ t <-*V 



e 2 = 



\/2tto- 2 



By (2.22) and the self-similarity of (Z t ) t > , f° r V > 0, 

1 - e-^ y 



/ i i \ 1 — e 



-Viy 



V 



2G 



which, when multiplied by exp(— y 2 /2cr 2 ), becomes integrable on [0, oo). Hence, by (2.21) and the dominated conver- 
gence theorem, 
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(B.17) 



To find the asymptotic behavior of the first integral in ( |B.16 ), we decompose it as: 



E (e- u *-l)l. 
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u t 



e~ u du 



(B.18) 



For J2i(t,y), using Markov's inequality and (3.2 ii), there exist k < oo and < to < 1 such that 

J2i(t,y) = f e- x v(z t < -tiy,t-YU t < t~^xjdx 



<E(e"' 1/YUt ) I e~ x exp (t~^x)dx < k- 

J — oo 1 



t.y 



for any < t < to and y > 0. Since 1 — Y/2 < 1/2 < 1/Y, for 1 < Y < 2, by the dominated convergence theorem, 
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(B.19) 



We further decompose the second term J22(t,y) in (B.18) as: 

J22(t,y) = E((e' Ut - 1 + U t )l t i ~ A-E(u t l, 

K ' V- ' {Z t <-t2y,Ut>0}) V {Zt 

=:4$(t,y)-J$(t,y). 



<-t2y,U t >0} 



(B.20) 



Using (JJUji) and the fact that 1 - Y/2 < 1/2 < 1/Y, for 1< Y < 2, 

0<4f(t,y)<t^E(t-^ Y U t y 



= t^ Y 
<t 1 / Y (l + 



t- 1/Y U t > u)du 



P (t- l,Y U t > uj du 

- tl/Y ( i+ 1°° Rt ( tl/Yu y Ydu ) 



(B.21) 

(B.22) 

(B.23) 
(B.24) 



which is clearly o(t 1 2). Moreover, 
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dw. 



Using (B.10), (1 - e'^/Vt < y, for y > 0, and applying the dominated convergence theorem, we get 
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(B.25) 
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1-Y 



as t -> 0. Combining ( |B.17| , ( |B.19| , ( |B.24[ ) and ( |B.25[ ), lead to 

limt^hit) = — y 
Step 3. Wc finally study the behavior of Isit). Note that 



h(t) = e-V E( e *l r i 1 / 



(B.26) 
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We first investigate J^\{t, y) := E I e Ut l 



{Z t >t2y} 

J 31 (t,y)=E((e-^-l)l 
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by decomposing it as: 



dy 
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{Z t >t2y} 

4 1 1 \t,y) + 4 2 1 ) (t,y). 



Z t > t*y 



(B.28) 



By (2.22), it is easy to see that J^(t,y) < nt 1 2 y Y , for any < t < 1 and y > 0. Hence, dominated convergence 



theorem together with (2.201 lead to 



lime-^ii- 1 / jg\t,y) 



1 - e^ y eri?- 



\ft V2TTCT 2 



dy 



9± 

Y 



^TTCT 2 



dy. 



(B.29) 



As usual, (t, y) can be further decomposed as: 



4[\t,y)=E((e- u '-1)1 



'{Z t >t3y,C/ t >0} 



Now, for u > 0, y > and £ > 0, by Markov's inequality and (3.2 -ii), there exist k < 00 and < to < 1, such that 
P^Zt > tiy, U t < ~uj < ¥(p t < -uj < E (e-^ 1 ^ ^ e-^ 1/Yu 



e u P(Z t > t2y,U t <-u)du 



< He 
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for all < t < to. Also, as we did in (B.24), there exist K < 00 and < to < 1 such that 



< E(7l - e- Ut )l r 1 ~ ,) <E(u t ) + <Kt 1/Y , 
for all < t < t . Hence, for any y > and < t < to < 1, and since 1 - Y/2 < 1/2 < 1/Y, for 1 < Y < 2, 



Jo 1 — *v 



(B.30) 



(B.31) 



(B.32) 



as t — > 0. Since both control functions in (B.30) and (B.31) are independent of y, combining (B.29) and (B.32), and 
by the dominated convergence theorem, 
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Next, by the self-similarity of (Z t )t>o, we decompose the expression Jm(t,y) :=E(e Ut l i (1 — e Zt )j\/t 
appearing in (B.27I as: 



J 3 2(t,y)=E[ 1 
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,-v t _ p-{Z t +U t ) 
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(B.34) 



(2") 

Note that J^ 2 , ( t, y) is quite similar to Ji2(t, y) in (B.4) and, thus, the corresponding integral has an asymptotic behavior 
similar to (B.6). Concretely, 
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Next, decompose J^l\t,y) as 
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Note that for x > 0, 
while for a; < 0, 
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Using the estimates (B.8| and (B.12), arguments as in getting (B.9) and (B.13I give 
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Combining (B.33I, (B.35 1 and (B.36I, 
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Finally, from ( |R2[ ), ( |B.15[ ), ( |B.26[ ) and ( |B.37[ ), and since 1 - Y/2 < 1/2, for 1 < y < 2, we obtain ( |42| ), therefore 
finishing the proof. □ 



Proof of Proposition |4.3[ 

When the diffusion component is nonzero, [121 Proposition 5] implies that ait) — > cr as f — > 0. In particular, att 1 / 2 — > 
as £ — » and, thus, (A.20) above remains true. Let cr(f) := at — cr, then <j(t) — > as t — > 0, and (A. 20 1 can be written 

as 
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By comparing (4.3)-(4.4) and (B.38), we have 

{C+ + C-)2- Y o 1 - Y 
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2tt 



t 1/2 , t^O, 



and, therefore, 
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The proof is now complete. 



□ 



C Additional proofs 



Verification of (2.22) 



It suffices to show that, for some constant k < oo and any v > and < t < 1, 

>^ i/y «) < K ti«r y . 

Indeed, using the same argument for —Zx, one can show an analog bound for P (Z\ < t~ x / Y v\ with v < 0, which in 



turn will imply (2.22). First, in light of (2.19), there exist < R < oo and H > 0, such that for any u> H, 
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for all < t < 1 and t> > 0. 



Verification of (B.5) 



This is similar to the proof of (2.22). Using (C.l), for any t > and y > 
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Then, using (2.22), 



+ t (^-D(i->')i ( RH i-y + H (_^L 

{t2-Yy<H} y \t2-7 



) 



) 



< r(y - i)-y- y + t^-h)^-y) Rt ih-^-Y) y i-Y 
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and (B.5) follows with A := RY/{Y - 1) + H Y . 



□ 
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